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^ • Abstract 

o 

00 ■ Jet-like three-particle azimuthal correlations can discriminate various physical sce- 
narios that have been proposed to explain the observed strong modification to two- 
particle azimuthal correlations. The three-particle correlation analysis is notoriously 
\^ difficult in heavy-ion collisions due to the large combinatorial backgrounds. We de- 
scribe the general idea behind the jet-like three-particle azimuthal correlation analy- 
sis, with emphasis on the subtraction of the combinatorial backgrounds. We discuss 
' in detail various sources of systematic uncertainties in such an analysis. 
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^ 1 Introduction 

Two-particle azimuthal correlations, between one particle (trigger particle) at high transverse mo- 
mentum (p±) to trigger on a jet and the other (associated particle) at lower p±, have become 
a powerful tool to study interactions between jets and the hot and dense medium created in rela- 
tivistic heavy-ion collisions at RHIC. The measured two-particle azimuthal correlations were found 
to be significantly modified by the medium [T1[2l[3|l]l5|6|7|l8|9l[T0] . This modification, in turn, can be 
used to study the properties of the medium. 

The away-side associated particles (opposite to the trigger particle in azimuth) in central Au+Au 
collisions, in particular, were found to be not narrowly distributed in azimuth at A0 = vr from 
the trigger particle, as observed in pp and d+Au collisions. Rather, the azimuthal distributions of 
the away-side associated particles are very broad, and perhaps exhibit a double-hump structure 
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with a dip at A0 = vr [2|3Pll5]|6]|71l8ll9][T0] . The away-side associated particles were also found to 
be significantly softened, not much harder than the bulk medium particles, indicating a partial 
thermalization between jet particles and the bulk medium PI51IHI51I71I5] . Moreover, the away-side 
associated particles were found to be softer at = tt than in the double-hump regions PIMTIIS] . 



Several scenarios have been proposed to explain the observations. One scenario is defiection of 
jets by transverse radial fiow of the bulk medium [11] or by the large survival probability of jet 
particles moving outwards than inwards due to the path- length dependence of energy loss |12j . 
Such a scenario would have jet particles narrowly clustered in individual events but the cluster 
is broadly distributed around A0 = tt over many events. Another scenario is large angle gluon 
radiation [13]. Such a scenario would possess quahtatively similar structure as for defiected jets. A 
third scenario is conical fiow generated by sound shock-waves due to large energy deposition by high 
momentum partons in the medium [13|15lll6j . Such shock- waves result in a distinctive Mach-cone 
type structure where particles are preferentially emitted at the Mach angle determined by the speed 
of sound in the medium, independent of the particle p±. If Mach-cone type conical fiow is indeed 
responsible for the observation, then the extraction of the speed of sound (thereby the equation 
of state) of the created medium may be possible. A fourth scenario is Cherenkov gluon radiation 
generated by interactions of fast particles with the medium [T71IT8] . Such a scenario would also result 
in conical emission of particles, as same as Mach-cone shock-waves, however, the Cherenkov angle 
would depend on the associated particle p± [IB] . 

Two-particle correlations cannot distinguish these scenarios because they give, quahtatively, the 
same two-particle correlation function. Three-particle jet-like azimuthal correlations are needed to 
distinguish these scenarios by using the extra information offered by the third particle. If the broad 
two-particle correlation is due to defiected jets or large angle gluon radiation, the two associated 
particles will be narrowly clustered in angle but the cluster will swing over a wide range in azimuth 
on the away-side. If Mach-cone shock-waves or Cherenkov gluon radiation is responsible for the broad 
two-particle correlation, then particles will be emitted conically, thus the two associated particles 
will have equal probabilities being opposite away from A0 = tt and being clustered together. A three- 
particle jet-correlation signal with opposite azimuthal angles from A(f) = vr for the two associated 
particles is therefore a distinctive signature of Mach-cone shock-waves or Cherenkov gluon radiation. 
The p± dependence of the cone angle will further discriminate between the two scenarios [TBJ. 

In central heavy-ion collisions, jet-like correlations are overwhelmed by large combinatorial back- 
grounds due to the large particle multiplicity. The signal to background ratio depends on the collision 
centrality and the p± range of the associated particles; the typical signal to background ratio for a 
modest associated p± range of 1-2 GeV/c is of the order of 1/30 for central Au-|-Au collisions. The 
three-particle correlation background will be much larger because the background level is squared; 
with the same associated p± range the three-particle signal to background ratio is of the order of 
10~^ in central Au-|-Au collisions. Thus the three-particle correlation analysis is extremely challeng- 
ing - one needs a statistical and systematic precision of the order of 10~^; the precise subtraction 
of background is critical. In this paper, we describe an analysis method of three-particle jet-like 
correlations [71J8] that is feasible to achieve such a precision. 

The paper is organized as follows. Section [2] describes the general ideas of the analysis method and 



2 



background subtraction technique. Section [3] discusses complications in a real data analysis and 
the handling of backgrounds. Section |4] discusses the various sources of systematic uncertainties 
involved in such a data analysis and demonstrates that the systematics are well controlled with the 
background subtraction technique. Section summarizes the paper. 



2 Analysis Method 



We suppose an event is composed, besides the trigger particle, of two components, the signal that is 
correlated with the trigger (such as a di-jet that is selected by the trigger particle but excluding the 
trigger particle itself) and the underlying background from the bulk medium (anything uncorrelated 
to the trigger particle except correlations due to anisotropic flow). Suppose the number of jet- 
correlated particles is A^assoc, the number of uncorrelated background particles is A^bkgd, and the 
total number of particles is = A'^assoc + A^bkgd- The total number of particle pairs is A^(A^ — 1), 
which is composed of three parts: the number of background pairs iVbkgd(A^bkgd — 1), the number 
of jet-correlated pairs A^assoc(A'assoc — 1), and the number of jet particle and background particle 
pairs (cross pairs) Ai'assocA^bkgd + A^bkgdA^assoc = SAi'assocA^bkgd- To obtain the number of jet-correlated 
pairs, one may subtract the numbers of background pairs and cross pairs from the total number 
of pairs [!]• Of course, if one is only interested in the number of jet-correlated pairs, one can easily 
get that from the number of jet-correlated particles, N^^soc- Since we are interested in the angular 
correlation of the jet-correlated pairs that is unknown a priori, the only way to do is to subtract the 
angular correlations of the background pairs and the cross pairs. The angular correlations of the 
background pairs and the cross pairs can be obtained because they, unlike those of jet-correlated 
pairs, are not intrinsically correlated (except for the anisotropic flow correlation). Therefore, by 
subtracting from the raw three-particle correlation the angular correlations of the background pairs 
and the cross pairs, one obtains the genuine three-particle jet-correlation. 

Below we are restating the above with a set of equations. The two-particle jet-correlation is 



y2(A0) = 1^2(A0)-E2(A0); (1) 
Y2iA(P) = dN/dAcf>, (2) 
S2(A0) = rfiVbkgd/c^A(/., (3) 

y2(A0) = rfiVjet/c?A0, (4) 

where A0 = (p — (pt is the associated particle azimuthal angle (0) relative to that of the trigger 
particle {4>t)- We have used Y2, B2, and Y2 as shorthand notations for the raw particle azimuthal 
distribution (relative to the trigger particle azimuth) for the entire event (i.e., raw two-particle 
azimuthal correlation per trigger particle), the azimuthal distribution of the underlying background 



^ We have used fixed multiplicities of N^ssoc and A'bkgd (and N = A^assoc + A^bkgd) in the text, but the 
argument is also valid for average multiplicities except that the numbers of pairs are (A^assoc)^j (Abkgd)^) 
and {N)'^ for Poisson or Gaussian multiplicity distributions. 



3 



particles, and the azimuthal distribution of the associated jet particles (i.e., the two-particle jet- 
correlation per trigger particle), respectively. 



The three-particle jet-correlation that we point out is the correlation between the trigger particle and 
two associated jet particles. What we have is the raw three-particle azimuthal correlation function. 



r3(A0i, A02) = l"2(A0i) ® r2(A02) = ^2(A0i) + 52(A0i) r2(A02) + 52(A0 



(5) 



where A0i = (f)i — (pt and A02 = <f)2 — 4't are the azimuthal angles of the two associated particles 
(01 and 02, respectively) relative to that of the trigger particle. We have used the notation of A®B 
to indicate correlation between A and B, not a simple product. The raw three-particle correlation 
function contains three pieces, the genuine three-particle jet-correlation that we discuss, 

%{A<f)u A02) = F2(A</)i) ® F2(A02), (6) 



and two background terms. One is the combinatorial background between an associated jet particle 
and a background particle, 

Y2 0B2 = F2(A0i) ® 52(A02) + 52(A0i) ® r2(A02), (7) 



which we term as the /iard-so/t background. If the jet-correlated particle and the background particle 
are uncorrelated, then we can remove the symbol from Eq.([7]), and readily obtain the hard- 
soft background from the simple product of the two-particle jet-correlation function Y2 and the 
background B2 as 

Y2 0B2 = r2(A0i)52(A02) + 52(A0i)F2(A02). (8) 



However, induced correlations between jet-correlated particles and background particles can be 
expected, for instance, due to their correlations both to the orientation of the trigger particle relative 
to the reaction plane angle and to the event centrality (impact parameter). Thus in principle the 
hard-soft background can only be obtained by 

y2(A0i; 0t, ^, b, ...)i?2(A02; 0t, ^, b, ...)+ 
fi2(A0i; 0t, ^, ...)^2(A02; 0t, ^, b, ...) 



00 2iT 2n 

Y B - [ db I f 

00 



where da and atot are the differential and total cross-sections, respectively, b is the impact param- 
eter, and ^ is the reaction plane angle (impact parameter direction). There could be other event 
characteristic variables besides b and \E'. The handling and effects of those induced correlations will 
be discussed later in Secions 14.21 and 14.31 respectively. For the sake of easy description and rela- 
tively simple equations, we will assume, for now, no correlation between jet-correlated particles and 
background particles, and carry on our study using Eq.([H]). 
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The other background term is the combinatorial background between two background particles, 



53(A0i, A02) = 52(A</)i) ® 52(A02), (10) 

which we term as the soft-soft background. The soft-soft background can be obtained from inclusive 
events as we will discuss in Section I3.3[ Once the two background terms are determined, the final 
three-particle correlation is simply obtained by 

r3(A0i, A</.2) = r3(A0i, A02) - [F2(A0i)52(A02) + S2(A0i)y2(A02)] - fi3(A0l, A02). (11) 
2.1 Elliptic Flow Correction 



In heavy-ion collisions, the background particles are correlated to the reaction plane due to the 
hydrodynamic kind of collective flow of the bulk medium and the anisotropic overlap region between 
the colliding nuclei. The trigger particle emission is also correlated to the reaction plane due to, not 
so much of the hydrodynamic collective flow, but the path-length dependent energy loss of high p± 
particles in the medium that is initially anisotropic. This reaction plane correlation at mid-rapidity, 
expressed in harmonics up to the fourth order, is given by 

= ^[1 + 2v, cos 2(0 -^) + 2v, cos 4(0 - v]/)], (12) 

where is the azimuthal angle of the particle and \E' that of the reaction plane. The reaction 
plane distribution is uniform. The flrst and the third harmonic terms are, due to symmetry, zero at 
mid-rapidity for symmetric collision systems which we consider in the present work. 

The distribution of the number of trigger-background particle pairs is given by 



d^N _ 1 Ntrig 



d'^d(f)td(f)dA(j) 271 2n 
' 2tt 



1 + 2ff cos2(0t - ^) + 4t;J^ cos4(0t - *) 
[1 + 2^2 cos 2(0 - ^) + 4^4 cos 4(0 - ^)] • 5( A0 - + 0t) (13) 



for a given set of the reaction plane angle the trigger azimuthal angle 0t, and the associated 
particle azimuthal angle 0. Here A^trig is the total number of trigger particles, A^bkgd is the number 
of background particles per event, f 2 and V4 are the second (elliptic flow) and fourth harmonic coef- 
flcients of the background particles, and the superscript '(t)' denotes those of the trigger particles. 
Integrating over 0t, 0, we obtain 

dN 'r,/r,,'r,, d^N 



rfA0 J I I rf^rf0trf0rfA0 
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-/VtrigA^bksd 



27r 



1 + 24*^172 cos(2A0) + 2v^4^Vi cos(4A0) . 
Normalized per trigger particle, the anisotropic flow modulation is therefore 
B2{A(j)) = Si [l + 24*^^2 cos(2A(/)) + 2vf\4 cos(4A0)] , 

where Bi = A^bkgd/27r is the background azimuthal density. 

The distribution of trigger-background-background triplets is given by 



1 iVtrig 



d'i/d4>td4>id4>2dA4>idA(l)2 2n 2n 

' 2n 

N2 

' 2% 

■ S{A(l)i - 01 + (/.t) • S{A(I)2 - 02 + 0t) 



1 + 2v^2^ cos 2(0t - + 4vf ^ cos 4(0t - 
1 + 2vi^^ cos 2(01 - + 4vf ^ cos 4(01 - 
1 + 2^^^ cos 2(02 - ^) + 4vf ^ cos 4(02 - ^) 



for a given set of 0t, 0i, and 02- Integrating over 0t, 0i, and 02, we obtain 



(14) 



(15) 



(16) 



2lT 



2tt 



2tt 



2tt 



■ j d"^ J d(j)tj d<pi j d(j)2 



d^N 



NtrMN2 



X 



dA0idA02 J 'J '"^J '"'I d*d0trf0id02dA0idA02 (27r)2 

^ 1 +2v^2^vi^^ cos(2A0i) + 24*^4^^ cos(2A02) + 2v^^Ki^^ cos2(A0i - A02) ^ 
+2^1*^1;^^ cos(4A0i) + 2v^'^vf'^ cos(4A02) + 2i;f ^i;f ^ cos4(A0i - A02 



+2v^2^vi^^vf'^ cos2(A0i - 2A02) + 2v^2^ v^2^ v^-^ cos2(2A0i - A02) 

\^ +24^^^^^^ cos 2(A0i + A02) 



;i7) 



Here we have taken the two background particles to be from two separate sets, with A^i and N2 the 
respective multiplicities. The product N1N2 is the number of pairs. For particles from the same set, 
this becomes (A^bkgd(-^bkgd ~ !))• Normalized per trigger particle, the anisotropic flow modulation 
is therefore, 

i^3(A0i,A02) = i^^^^i^5^^ 

( 1 +2vf cos(2A0i) + 2wf cos(2A02) + 2vi^\!^'^ cos2(A0i - A02) ^ 
+2^1*^ wi^^ cos(4A0i) + 2t;l*^wf ^ cos(4A02) + 2v^4'^vf cos4(A0i - A02) 

(18) 

+2v^\!^^vP cos2(A0i - 2A02) + 24*^4^^-^?^ cos2(2A0i - A02) 
\^ +2^'^^;^ cos2(A0i + A02) j 
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2.2 Hard-Soft and Soft-Soft Backgrounds 



The two-particle jet correlation signal can be obtained from Eq. ([T]) and ([T^ as 

f2(A0) = Y2iA(f)) -Bi\l + 2v^2^V2 cos(2A(/)) + 2v^\i cos(4A(/))' . 



(19) 



The hard-soft background is given, as in Eq.®, by folding the two-particle jet-correlation signal of 
Eq. flTIJl) with the underlying background of Eq. flT^ . namely 



1^2 ® ^2 = iY2{A(f)i) -Bi 1 + 2v!2^vi^^ cos(2A0i) + 2v^^>vi'' cos(4A0 



Bi [1 + 24*^ff ' cos(2A02) + 2vf'vl^' cos(4A02 
(y2(A02) - 5i [l + 24*)t;f cos(2A02) + 2t;f cos(4A(/.2)]) 
(5i [l + 24*^4^^ cos(2A0i) + 2ff ^ cos(4A0i) 



(t)„,(2) 



+ 



X 



X 



(20) 



The soft-soft background is given by Eq. f|T8|) . Obviously the anisotropic flow is the only considered 
correlation in the soft-soft term of Eq. f|T8l) . There are other correlation effects in the real collision 
data, which we will discuss below in section 13.31 



3 Complications in Real Data Analysis 



In real data analysis, the events are often binned in centrality. The centrality is often defined in 
terms of particle multiplicity. The background shape can be obtained from mixed-events, where 
the trigger particle is mixed with inclusive events LH from the same centrality bin; it can be also 
obtained from single particle distributions of minimum-bias events folded with a random trigger 
particle. These two should give the same result, however, the former may be preferred if subtle two- 
particle acceptance effects are present that cannot be taken care of in the single particle folding. 

Before event-mixing or particle folding, the particle azimuthal acceptance has to be corrected such 
that the azimuthal distribution average over many events is uniform. This correction should be 
applied to both trigger particles and associated particles. It is equivalent, in the case of two-particle 
correlation analysis, to the two-particle acceptance correction as is usually done in data analysis. 
For three-particle correlation, however, this cannot be simply corrected for by the folding of two 
two-particle acceptance corrections. 

The analysis of three-particle correlation is extremely challenging because of the large combinatorial 
backgrounds. For central Au-|-Au collisions, the number of charged particles at mid-rapidity within 

^ The term "inclusive events" refers to events from the minimum bias event sample within the same 
centrality class as those with a trigger particle but without specifically requiring a trigger in them. 
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\rj\ < 1 for a modest j9_l range of 1 < < 2 GeV/c is about 150; the di-jet charged hadron 
multiphcity within the same range for a trigger particle range of 3 < < 4 GeV/c is of 
the order of a few (~ 5) [7IITU] . In other words, the signal to background ratio is of the order of 
1/30. Thus, the hard-soft background is a factor of 30 larger than the three-particle correlation 
strength, and the soft-soft background is a factor of 10^ larger. In order to have a significant three- 
particle correlation result, one needs a precision of lO"'', both statistically and systematically. Such 
a precision is unprecedented in relativistic heavy-ion collisions; the background subtraction is crucial 
in three-particle correlation analysis. 

3.1 Two-Particle Correlation Background 

In real data from relativistic heavy-ion collisions, because exclusive jet reconstruction is impossible at 
low jet energy, the background level and the true jet signal strength are unknown a priori. One has to 
make an ad hoc assumption about the background level and hence the signal strength. The often used 
assumption is that the jet-correlation signal is zero at some A0 region. The two-particle azimuthal 
correlation function measured at RHIC is the lowest around = 1. The STAR experiment makes 
the assumption that the jet signal is zero within the fixed range of 0.8 < |A(/)| < 1.2 (which we will 
refer to as zero- yield- at- 1, or ZYAl) [2]. The PHENIX experiment uses the so-called zero-yield-at- 
minimum (ZYAM) method in which the A0 region where the signal minimum resides is determined 
by the data itself [5j. 

There are many sources of background, including the hydrodynamic type underlying event consisting 
of mostly soft particles, as well as hard particles from other jets not selected by the trigger particle. 
Since the underlying background event is uncorrected with the di-jet selected by the trigger particle, 
the underlying background event should be represented by the inclusive events (with a caveat 
discussed below). This is not obvious at the first glance because, for example, the underlying event 
can contain other jets that are not selected by the trigger particle as aforementioned. But so do 
the inclusive events. We can assure ourselves by realizing the following. Suppose the number of jets 
(with a trigger particle in the selected p±_ range) is Poisson distributed with an average of (n) in 
inclusive events, the probability to have n jets per event is Pn = (n)"e~^"'^/n!. We select events with 
at least one jet (with a trigger particle in the selected p± range). The probability to have n + 1 
such jets (that also contain trigger particles) in the same event is Pn+i = {n)'^^^e~^'^'^ / (n + 1)!. 
This is equivalently the probability to have a triggered event (with a trigger particle) with n 
jets (that have trigger particles) in the underlying background event that is not associated with the 
selected trigger particle. But this event is going to be used n + 1 times (because all trigger particles, 
not just the highest p± particle, are used in the analysis), so the probability to have n jets in the 
underlying background event, now normalized by the total number of trigger particles, is given by 
{n + l)Pn+i/ Z]^o('^ + ^)Pn+i = {n)''^e~^^'^ /n\. This is identical to P„, the probability to have n jets 
per event for the inclusive events. 

Since the inclusive events represent the underlying background, the background shape can be con- 

^ We use "triggered event" to refer to the event selected by a high pj_ trigger particle in off-line data 
analysis. We do not mean an event triggered on-line. 
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structed by mixed-event technique, mixing trigger particles with all particles (within the associated 
particle p± range) from the inclusive events of the same centrality. Because the event plane corre- 
lation is lost in such event-mixing [20j, the flow modulation needs to be added in by hand. 

There is one caveat, though, in using the inclusive events to represent the underlying background 
in triggered event. This is often referred to as "trigger bias": because of the centrality cuts often 
used in data analysis, and because the production of soft particles is proportional to the number of 
participants whilst that of high-p^ particles is proportional to the number of binary collisions (but 
is strongly suppressed in central heavy- ion collisions), the inclusive events (that are from the same 
centrality class as those triggered events) do not exactly represent the multiplicity of the underlying 
background event. Thus, the background level needs to be scaled with the aforementioned working 
assumption of two-particle correlation ZYAl (or ZYAM): 

Bi = aB,^,. (21) 



Here the scaling factor a is the ratio of the underlying background multiplicity (in the triggered 
events) to the inclusive event multiplicity for a given centrality: 

a = A^bkgd/A^inc- (22) 



The background estimated from the two-particle ZYAl (or ZYAM) method is likely an upper limit 
because the jet signal, due to its broadening, may be positive in all azimuth. As will be discussed 
in Section 14.51 the three-particle correlation contains richer information because of the presence 
of a third particle. One may require the three-particle correlation to be positive definite and use 
three-particle correlation ZYAM to obtain the scaling factor a. The obtained a will be smaller than 
that obtained from two-particle ZYAM, and may be more accurate in representing the underlying 
jet-correlation background. 



3.2 Hard-Soft Background in three-particle Correlation 



As discussed above we use aBi^c to estimate the multiplicity level Bi of the underlying background 
in triggered events. The two-particle correlation is thus given by Eq. (HM with Bi replaced by aB^^c'- 



r2(A0) = F2(A0) - aAnc 1 + 2vi'\2 cos(2A0) + 2v)^'v4 cos(4A0) 



(23) 



Likewise the hard-soft background is given by Eq. fl20l) with Bi replaced by aB^ 



Yo^Bo 



aBu 



1 + 24*^4^^ cos(2A0i) + 2<'<^ cos(4A0 



(t).,(i) 



(aBi 
Y2{A<f)2 



1 + 24*^4^^ cos(2A02) + 2<^<^ cos(4A02 



(t).,(2) 



+ 



aBu 



1 + 24*^4^^ cos(2A02) + 2<^<^ cos(4A02 



(t).,(2) 



(aAnc [1 + 2f f cos(2A0i) + 2vf\^l^ cos(4A0i 



(24) 
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In real data analysis, the background is often obtained from mixed-events, mixing trigger particles 
with inclusive events as mentioned above. The event-mixing technique is used to account for any 
detector effects, which are present in the correlation signal. The flow modulation is often added 
dynamically pair-by-pair to avoid potential systematic effect of {v2^V2) ^ ('f^2*^)('f^2)- However, such 
a systematic effect is found to be very small in real data. 



3.3 Soft-Soft Background in three-particle Correlation 



The soft-soft background can be obtained by mixing inclusive events with a random trigger (from 
a triggered event), expressed in two (p angles relative to the random trigger particle: 

B, (A0„ A02) = ^ J d^, ^^^^A^^^A^^ • (25) 



This is effectively the two-particle distribution 02) in the two azimuthal angles but with the 

origin of the azimuthal angles shifted to the random trigger particle direction. This event-mixing 
procedure takes into account subtle acceptance effects that may not be fully corrected for by the 
single particle acceptance corrections. If the acceptance corrections are perfect so that the corrected 
single particle distributions of trigger and associated particles are perfectly uniform in azimuth, 
then the event-mixing is unnecessary and i?3°'^(A0i, A02) = i?3°'^(0i, ^2). 

In real events, the background particles contain not only anisotropic flow correlation, but also other 
correlations that are unrelated to the trigger particle, such as jet-correlations due to jets other than 
the di-jet selected by the trigger particle [21] ■ Those correlations (other than the flow correlation) 
are contained in the inclusive events used in the event-mixing, therefore are already included in the 
soft-soft background in Eq. (125|) . 

The anisotropic flow correlation between the two soft particles are also already contained in the 
inclusive events being mixed, but the flow correlation between the soft particles and the trigger 
particle is lost because of the randomization of the reaction planes in the mixing [2U]. The trigger- 
flow correlation needs to be added in by hand for each trigger-background-backgrounbd triplet in 
the event-mixing via Eq. f|T8|) except the two terms that do not contain the trigger particle flow. 
The trigger-flow modulation that needs to be added in by hand is given by: 



B 



X 



i,,,tf _ (iVmc(iVinc - 1)) 

(2vr)2 

2t;f cos(2A0i) + 24*^ff ^ cos(2A02) + 2v!i\^^^ cos(4A0i) + 2vPv!i'^ cos(4A02) + 
2t;f ^ cos2(A0i - 2A02) + 2ff ^ cos2(2A0i - A02) + 
y 2vi^K^2^vf'^ cos 2(A0i + A02) 



(26) 
\ 

/ 
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The trigger-flow correlation can be constructed by mixing the trigger particle with two inclusive, 
within the same centrality window but otherwise random, events. 

The multiplicity of the inclusive events, (A^mc)? does not exactly represent the underlying background 
multiplicity of the triggered events, (A^bkgd)- Thus, we scale the soft-soft background constructed 
from the inclusive events, given by the sum of Eqs. fl2^ and (1271) . by the ratio of the number of 
pairs in the underlying background to that in the mixed inclusive event: 

(A^bkgd(A^bkgd - 1)) _ 2 (cyrj. 



where a ~ 1 is given by Eq. (!22|) . Eq. (!27|) would be exact if the multiplicity distributions are 
Poisson (or Gaussian in the large multiplicity limit) for which {N{N — 1)) = (A^)^. Even if the 
multiphcity distributions are not Poisson, as long as ^^"""^^^^^p"^^^ = ^^'"f^^a"^^^ ^ Eq.(|27D would 
still be exact. However, in real data, the multiplicity distributions in the inclusive events and the 
underlying background of the triggered events can deviate differently from Poisson distributions. 
To take this into account, a second scaling factor b is applied in front of a^, namely 



(A'bkgd(A^bkgd — 1)) _ j,^2 



(iVinc(iVinc - 1)) 



ba\ {2t 



Here the factor b characterizes the effective difference in the associated multiplicity distributions of 
the inclusive events and the underlying background of triggered events: 

, ^ (iVbkgd(iVbkgd - 1)) /(iVbkgd)^ . . 

(iVinc(iVinc - 1)) /(iVinc)^ ' ^ ^ 



The value of b should be very close to 1. Due to the large combinatorics, even a small change in b 
could have a significant effect on the magnitude of the final three-particle correlation result. However, 
the shape of the three-particle correlation is not much effected by small changes in b because the soft- 
soft background is largely uniform. Since triggered events in heavy-ion collisions are dominated by 
the underlying background, we assume the characteristics of the associated multiplicity distributions 
in the entire triggered event and in the underlying background to be the same, ^^bkgd(A^bkgd- i)) _ 

^'^^^jT^, and obtain b as 



{NiN-l))/{Ny 

(iVinc(iVinc - 1)) /(iVinc)^ 



(30) 



Thus the complete soft-soft background is given by 



5;, = ba' 



5r(A0i, A02) + bT'\A(P,, A02) 



(31) 
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The trigger- flow by Eg. (1271) would be exact if the background particles were only correlated by 
flow. However, additional correlations are present between the background particles in real data 
as already included in Eq. fl2H) . such as those due to resonance decays and background (mini-)jets. 
Those additional correlations alter the trigger-flow background. The effect can be corrected for to 
the flrst order by applying a multiplication factor Bi^'^ / B^^'^'^^ in front of B^^'^'^^ given by Eq.(I2j 



Here 5™'^'^^ would be the soft-soft background in analogy to B™'^ of Eq. fl25l) if all correlations were 
due to flow. It can be obtained from the measured soft particle flow, keeping the number of soft 
pairs as same as in -B™'^, by 



(iVinc(iVinc - 1)) 

(27r)2 



"l + 2vi^^vi^^ cos2(A0i - A02) + 2v^^\^^^ cos4(A0i - A02^ 
Therefore, the complete soft-soft background in Eq. (l3T!) is modifled as 

5r*'(A0i,A02r 



Bs = ba'Bl^'{A(f)i,A(j)2 



1 + 



B'r\A<p„A<p,) 



(32) 



(33) 



And flnally the three-particle correlation function is given by Eqs. (ITTI) . ( l24l) . and ([33 



4 Sources of Systematics 



In this section, we discuss the various sources of systematic uncertainties involved in the three- 
particle jet-like azimuthal correlation analysis. This is not meant to be a complete list of systematics, 
but we think these are the major ones. We conclude that the dominant systematic uncertainties come 
from uncertainties in the elliptic flow measurements and background normalization as discussed in 
sections 14.11 and 14.51 respectively. 



4-1 Elliptic Flow Uncertainties 



The anisotropic flow correction has effect on the flnal three-particle azimuthal correlation function, 
because the number of extracted associated particles depends on the amount of flow subtraction [2]. 
The size of the effect on the three-particle correlation function is similar to that on the two-particle 
correlation function. The sizes of the flow corrections in the two background terms (hard-soft and 
soft-soft) separately are, however, large compared to the magnitude of the genuine three-particle 
correlation, but the leading terms of the flow corrections in the hard-soft and soft-soft backgrounds 
cancel. We illustrate this by examining the total background, sum of Eq. ( l20l) and ([TF 



Y2 0B2 + B3 = F2(A0i) X 5i 1 + 2v^2 



•^*^f5^^ cos(2A02) + 2ff ^ cos(4A02) 



+ 
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r2(A(/)2) X Bi[l + 24*^4^^ cos(2A0i) + 2v^^\^^^ 



2Bl 



1 + 2V^2^V^2^ COS 



;(2A</)i) + 2vfv^l^ cos(4A0i)' 
^^^'-y^-' cos(2A(/)2) + 2vfvf^ r^.(AAA.\ 



I 



cos(4A02) 

1 +24*^4^^ cos(2A0i) + 2ff ^ff ^ cos(2A02) 
+2vfv^^ 

+ 2V2^V^\f^ COS 



cos(4A0i) 

X 

+ 



^X'vX" cos(4A0i) + 2vfvf^ cos(4A02) 



V + 



2(A0i - 2A02) + 2v^2^vf\^l^ cos2(2A0i - A02 
-2v^2^vf^vf cos2(A0i + A02) 



(34) 



Here we have dropped the anisotropic flow correlation between the two soft particles because they 
are already included in the inclusive events used in constructing the soft-soft background. In other 
words, the flow between the soft particles has been precisely included by the background inclusive 
event itself. We have also taken the number of background pairs per (2-7?)^ to be B\ for simple 
illustration. Rearranging the terms gives 



Bs = Y2{A<f)i)xBi 

F2(A</.2) xBi \l + 2v^^v^^ 



1 + 2i;^*^4^' cos(2A02) + 2v^^\P 



':2'v^^' cos(2A0i) + 2vt\^4^ 



I 



Bl 



+2v^2^v^2^ cos(2A0i) + 2ff cos(2A02) 
+2vfv^^ 



jX'v''^' cos(4A02) 
)X'vX' cos(4A0i) 



+ 



"4 V, 

-2v'i\^2^vf^ cos2(A0i - 2A02) - 2v'i\^2^ v^l'^ cos2(2A0i - A^a) 
-2v'i\f^vf cos2(A0i + A02) 
+S{vffv'i\P cos(2A0i) cos(2A(/)2) + Svfvfv'i-\f^ cos(2A(/)i) cos(4A(/)2) 
+S{vffv^l\f^ cos(4A0i) cos(4A02) + Svfvfvf\^l^ cos(4A0i) cos(2A02) j 
= 5i[r2(A0i) + F2(A02)-5i] + 
5i[r2(A(/)i) -5i]2t;f^t;^^^cos(2/ ' 



cos(4A0i) + 2vfvf^ cos(4A02) 



V 



(2v^\^^ 



Bl 



i^i[r2(A0i) + F2(A02)-i^iJ + 

5i[r2(A(/)i) - Bi]2v'i\f'^ cos(2A02) + 5i[r2(A(/)2) - Bi]2v'i\'i^ cos(2A(/)i) 
5i[r2(A0i) - Bi]2vfvf^ cos(4A02) + 5i[r2(A02) - Bi]2vfv^l^ 

j^2^v'i\f^ cos2(A0i - 2A02) + 2v^2^v'i\^l^ cos2(2A0i - A02)+ 



'4 cos(4A(/)i) 



^?v^2^vf cos2(A0i + A0; 



i(2A0i) cos(2A02) - Svfvfv'i^f'^ cos 



U2 V2 cos I 

U4 (^4 COS 



;(4A0i) cos(4A02) - Sv'i\fvf\^l^ cos 



i(2A0i)cos(4A02)- 
i(4A0i) cos(2A02) 



(35) 



As can be seen, the leading term in the background is the first term (first line of Eq. (135|) r.h.s.), 
of the order B^. The terms on the second line of Eq. (135|) r.h.s. are of the order l^i^i^g. Since Y2 
is of the order a couple of percent of Bi [71IT9] and is of the same order [22j, those terms are of 
the order or Y2 ~ Is. The terms on the third line are of the order -B^ff because ~ | 23j . 
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and those on the fourth hne are a mixer of Bfv2, B1V2, and smaller terms. Neglecting the terms of 
the order Bfv2 or smaller, we obtain 



Y2®B2 + B3 = Bi[Y2{A<f)i) + r2(A02) - 5i] + 

2Bi[Y2{A^i) - cos(2A02) + 2Bi[Y2{A(f)2) - Bi]vi'^vi'^ cos(2A0i) + 

2Bfv!2^vi^\f^ cos2(A0i - 2A02) + 2Blvi'K^2^ v^^^ cos2(2A0i - A02) + 
2Bfv^^\i^\^2^ cos2(A0i + A02) - 85^(4*^)^^^^^;?^ cos(2A0i) cos(2A02). (36) 

Except the first, leading term, all other terms in Eq. fl36|) are of the same order of the three-particle 
correlation strength, ^3. Therefore, uncertainties on V2 will result in uncertainty on the three-particle 
strength of the similar magnitude. 

In the above illustration, we have taken Bi to be constant, not affected by the uncertainties on V2- 
However, with the background normalization scheme of ZYAM or ZYAl, Bi will be affected by the 
value of ^2 P]- This effect on Bi, being obtained by normalization fixed at A0 = 1, works in the 
opposite direction as the V2 uncertainty effect on the three-particle correlation strength, thereby 
making the effect of V2 uncertainty smaller than illustrated above. 

Experimentally, anisotropic flow is measured within flnite centrality bins. The flow magnitudes 
vary with event multiplicity within a given flnite centrality bin. Since flow correction is not linear 
in flow (for both two-particle and three-particle correlations), the calculated flow correction using 
the measured average magnitude of flow in a given centrality bin deviates from the true flow effect. 
The magnitude of this effect can be estimated by parameterizing the measured flow as a function of 
event multiplicity. With the measured multiplicity dependence of flow, the effect is very small and 
can be neglected [21]. 

The major uncertainties on flow measurements are due to non-flow [22]. The tranditional reac- 
tion plane measurements have large contributions from non-flow, so do the two-particle cumulant 
measurements. The non-flow effect is signiflcantly reduced in the four-particle cumulant measure- 
ments because the jet multiplicity in our relavent jet energy range is small, mostly fewer than four 
particles. On the other hand, flow fluctuations reduce the measured four-particle cumulant V2 and 
increase the reaction plane and two-particle cumulant V2 [2S]- -Due to the presence of other jets 
not correlated with the di-jet selected by the trigger particle, non-flow effect should be partially 
included in the flow correction in jet- correlation measurements. So the V2 used in the correlation 
background should include the flow fluctuation effect but should not include the full non-flow effect. 
The interlink between flow and jet-correlation is complex; in the lack of any flrmer solution, the 
average of the four-particle cumulant V2 and the reaction plane V2 has been used for background 
correction in jet-correlation measurements [2]. The systematic uncertainty boundaries are set by 
the four-particle cumulant V2 and the two-particle cumulant V2- Here the two-particle, not the re- 
action plane V2 is used as the systematic limit because the flow fluctuation (which is present in 
two-particle V2) should be included in the background correction. However, in practice these two V2 
measurements are approximately equal. 
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4-2 Systematics due to jet anisotropy 



The hard-soft background, given by Eq. fl2Ul) or fl2^ . arises from the combinatorial pairing of a 
jet-correlated particle and a background particle. The background is modulated by anisotropic flow 
because the trigger particle and the background particles are all correlated to the reaction plane. 
In Eqs. ( 120|) and (^^, the measured trigger t>2 is used. However, the anisotropy of the trigger 
particle, with the additional requirement that a jet-correlated associated particle present within 
the corresponding kinematic region, may not be equal to the measured anisotropy for inclusive 
trigger particles a priori. We will refer to those trigger particles with an associated partner as jet 
trigger particles, and refer to this effect as jet anisotropy. The other trigger particles (i.e., having no 
associated particles in the corresponding kinematic region) of course need to have anisotropy also 
different from the measured inclusive one, in order to have the inclusive one (including all trigger 
particles) as measured. 

Suppose the elliptic flow of the trigger particle accompanied by at least one associated particle is 
V2 , and the fraction of those trigger particles is a. Suppose the elliptic flow of all other trigger 
particles is ^2*°'*. The measured elliptic flow of inclusive trigger particles, \ requires: 

avt^ + {l-a)vt^ =v^l (37) 



The two-particle correlation normalized by the total number of trigger particles, given in Eq. f[T^ or 
(!23|) . is unchanged. The two-particle correlation for jet trigger particles, normalized per jet trigger, 
is however a factor of 1/a of that given by Eq. (fT9|) or (1231) . The two-particle correlation for the 
other trigger particles, of course, vanishes. Only the jet trigger particles contribute to the hard- 
soft background; the other trigger particles do not contribute because of the vanishing associated 
particle multiphcity. The hard-soft background is thus given by 



i (F2(A0i) - 5i 1 + 2t;f cos(2A0i) + 2ff ^f['^ cos(4A0i 



Yo^Bo 



a 



Bi 1 + 24*^^4^^ cos(2A02) + 2t;. 



(ta), ,(2) 



cos(4A02) 



+ 



I (F2(A02) - 5i 1 + 2vfvf^ cos(2A02) + 2<^<^ cos(4A02 



(t).,(2) 



(3^ 



Bi 1 + 24*^^4^^ cos(2A0i) + 2vl 



(ta)„,(l) 



cos(4A0i 



where the normalization is, as always, by the total number of trigger particles. The 1/a factor in 
the two-particle correlation strength per jet trigger particle is cancelled by the fraction of the jet 
trigger particles of a. Note that the ^2*'* ('^4*^) the background modulation (that is folded with 
the associated particle signal) in Eq. (l20l) and (p^ are replaced by f2*'*'* (f4*'^'') in Eq. (1381) . 

The soft-soft background, as given in Eqs. (133|) . fl27|) and (132|) . is unchanged, because the in 
Eq. fl27|) is replaced by 0^2*^-* + (1 — «)f2*°'* which is identical to ^2*'' by Eq. fl37|) . 

Thus, the only change in the three-particle correlation analysis, due to the possible effect of different 
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anisotropies of trigger particles with and without accompanying associated particles, is the hard- 
soft background term. The estimation of V2^^ is somewhat arbitrary. One may use a number of 
inputs to make an educated guess; these inputs include (i) two-particle correlation with respect 
to reaction plane, (ii) how large and small the elliptic flow of other trigger particles can have 
with considerations of hydrodjTiamic calculations, and (iii) how small the factor a can get without 
running into unreasonably large two-particle correlation strength per jet trigger particle. Another 
reasonable approach is to use, as an estimation of f 2 , the measured elliptic flow of particles at p± 
equal to the sum of those of the trigger particle and the accompanying associated particle. 

However, there are further complications. We know that the flow background varies with the trigger 
particle azimuthal angle relative to the reaction plane |26] . Only when the two-particle jet-correlation 
is constant, is the hard-soft background from folding as in Eqs. fl20l) and fl2^ valid. In other words, 
Eqs.f l20p and fl2^ have the implicit assumption of constant jet-correlation over the trigger particle 
orientation. Eq. flSSl) has loosened this requirement in the sense that the jet-correlation magnitude 
can change with the trigger particle orientation (but still regularly according to harmonic osillations) 
through the different flow harmonics, V2^^ 7^ ^2*''. In real data, however, the jet-correlation may 
change with the trigger particle orientation in both its magnitude and shape. Such a change cannot 
be simply described by a differing f g*'*^ value. The hard-soft background can only be constructed from 
folding of the jet-correlation signal and background at each trigger particle orientation relative to 
the reaction plane separately and then integrated together over all trigger orientations. This requires 
measurements of two-particle jet-correlation as a function of trigger particle azimuthal angle relative 
to the reaction plane |27j . 



4.3 Induced Correlation between Associated and Background Multiplicities 



The hard-soft background term is obtained by folding the jet signal {Y2) with the background {B2), 
both averaged over a given class (centrality) of events. The jet signal and the background level are 
not intrinsically correlated, but their magnitudes could be correlated via the event multiplicity. For 
instance, within a given centrality of events, if the jet signal increases with the event multiplicity 
(which in high multiplicity events is essentially as same as the background multiplicity), then the 
average of the product of the jet signal and the background level will not equal to the product of 
the individual averages of the jet signal and background level, i.e., {Y2B2) 7^ {Y2){B2)- We shall 
estimate how large the difference, 6 = ^ qq^j^ central heavy- ion collisions by using 

simple assumptions of jet multiplicity as a function of event multiplicity, N, within a centrality bin 
Ni<N < N2. 

If the jet multiplicity linearly increases with ^2 = 01 + C2{N/Niy/^, and since central heavy- 

ion event multiplicity is large, B2 ~ N/2tt — Y2 ^ N/2tt, the difference can be expressed as 5 = 
{N I )-{N){N I ) ^ _ c^/c2. Obviously if C2 = 0, i.e., there is no apparent correlation between 



Y2 and i?2, then 8 is zero. If we assume constant cross-section over N within Ni < N < N2 so 

1. X {6/7)(r7/3_i)/(^+i)_(3/4)(r4/3_l) 

we can easily calculate the various averages, then we have = — / , , ,/ — \ where 

^ ' (r-l)c+(3/4)(r4/3-l) 
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r = N2/N1. Let us plug in some reasonable numbers to get a feeling of how large the effect could 
be. If we take iVi = 400 and N2 = 500, then we get 6 = qISSo- ^1 = 0, then the relative 
difference is 5 ~ 0.14%. Given the two-particle correlation signal to background ratio hence the 
three-particle signal to hard-soft background ratio of 1/30, this difference is about 5% of the three- 
particle correlation signal. If ci = C2, then the difference is about a factor of 2 smaller. 

If the jet signal linearly increases with the event multiplicity, Y2 = Ci + C2N/N1, and again B2 ~ 
N/27r — ^2 ~ N/271', then the effect is larger because the apparent correlation between Y2 and B2 is 
stronger. The relative difference can be obtained as 5 = ^jN^^e]\^f^]v)) • ^S^i^ C2 = 0, the difference 
is zero. Again assuming constant cross-section over Ni < N < N2, we have 6 = 3(i^rX2^+i+r) • Take 
the same example, A^^i = 400 and A^2 = 500, we have 6 ~ "c+f 25 • ^'^^ ~ relative difference 

is 5 ~ 0.4%, which amounts to 12% of the three-particle correlation signal. If ci = C2, then again 
the effect is about a factor of 2 smaller. 

We should note that the above estimates are on the conservative side, because the away-side asso- 
ciated particle multiplicity has a weak dependence on the event multiplicity [2]. The dependence 
varies with the associated p±. At low p± the associated particle multiplicity increases with the 
event multiplicity; at high p± the associated particle multiplicity decreases; at intermediate p± the 
dependence can be very weak. 

For pp and d+Au collisions, due to the significant trigger bias, the correlation between the two- 
particle correlation strength and the background multiplicity via event multiplicity can be large. 
Our estimate indicates that the effect can be as large as 5 ~ 10% j24|. However, the hard-soft 
backgrounds in pp, d-|-Au, and peripheral Au+Au collisions are small, comparable to the three- 
particle correlation strength. So the effect is of the order of 10% of the three-particle correlation 
strength. One may reduce the effect by dividing the pp and d+Au events into several multiplicity 
bins. 



4-4 Centrality Bias of Triggered Events from Inclusive Distribution 



The underlying background in the triggered events may be different from the inclusive events be- 
cause the triggered events may be biased towards larger multiplicity due to the requirement of a 
high p\_ trigger particle. In other words, the multiplicity distribution of the underlying background 
differs from that of the inclusive events. This effect can be studied by first obtaining the soft-soft 
background from the inclusive events as a function of the event multiplicity in single increment, 
and then weighting those soft-soft backgrounds with the multiplicity distribution of the triggered 
events. The multiplicity distribution of the triggered events should be offset towards the lower end 
by the amount corresponding to the jet associated multiplicity; this offset can be done in various 
ways to assess the systematic effect. 

Our estimate indicates that the systematic effect due to the slight mismatch in the multiplicity 
distributions in central Au+Au collisions is small and generally negligible |24j|. The mismatch in 
the multiplicity distributions in the triggered events and inclusive (or minimum bias) events for pp, 
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d+Au and very peripheral Au+Au collisions can be large. However, since the soft-soft background 
in these collisions is substantially smaller, the systematic effects due to the multiplicity distribution 
mismatch are also not large [23]. Nevertheless, one may subdivide the minimum bias events from 
pp and d+Au coUisions in several multiplicity classes, to the extent that is statistically meaningful, 
to reduce this effect. 

Due to the different shapes of the multiplicity distributions, the flow correction obtained from 
inclusive events of a given centrality may also be slightly off. The anisotropic flow is measured 
in the predefined centrality windows, not for each event multiplicity value. Nevertheless, one may 
parameterize the measured flow as a function of the event multiplicity and assess the size of the 
effect it may cause. We estimate this effect to be very small and negligible |24j . 



4-5 Background Normalization 



In two-particle correlation analyses, one often normalizes the background to the minimum of the 
raw signal so that the signal at minimum is zero (ZYAM). The STAR collaboration normalizes the 
background to the raw signal in the fixed range of | A0± 1| < 0.2, about where the signal is minimal 
(ZYAl) [2]. The use of the fixed region is to avoid an additional running systematics, which for 
example could be affected by statistical precision of the data. One should note, however, this scheme 
of background normalization is only a working assumption, in lack of something better. Although 
the jet signal is expected to be non-negative everywhere in azimuth, the jet signal can be larger 
than zero in all azimuth if the jet peaks are broad and overlap, as likely the case in central Au-|-Au 
collisions. 



On the other hand, three-particle correlation with jet signal in the two-dimensional {A(j)i, ia<^2) 
angular space can provide a more stringent constraint on the background level, because the three- 
particle jet-like correlation should, in principle, have non-negative values anywhere on the two- 
dimensional (A01, A(/)2) plane. For example, when the three-particle correlation have some region 
where the signal vanishes, the projection onto one angular variable can in principle be above zero 
everywhere. In this case the background normalization to the minimum signal in the two-particle 
correlation over-subtracts the background, which could result in negative regions in the three- 
particle correlation function. 



Below we use a specific example to illustrate the above points. We define a jet-like two-particle 
correlation with a near-side peak and a broad double-hump away-side distribution: 



27ro"i 



exp 



2a! 



+ 



N2/2 



271^2 



exp 



(A0 + 



2al 



+ exp 



2af~ 



(39) 



We use 0"! = 0.4, a2 = 0.7, A^i = 0.7, N2 = 1.2, and 6 = 1. Obviously the correlation is positive 
everywhere. The minimum signal strength is at A0 1 and is about 0.12. It is not zero because 
of the overlap of the Gaussians. We suppose this jet-like two-particle correlation is presented in 
every event, and the background is fiat in azimuthal angle at Bq = dN/dAcj) = 150/27r. The raw 
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two-particle correlation is shown in the left panel of Fig. [H The solid horizontal line shows the real 
background level, and the dashed horizontal line shows the normalized background level such that 
the jet-correlation signal is zero at A0 ^ 1. Because of the overlap of the Gaussians, the normalized 
background is not the true background level but is overestimated. 




A(|) A(|)^ A(|)^ 



Fig. 1. Left panel: two-particle correlation signal sitting on a flat background (solid horizontal line). The 
dashed horizontal line is the overestimated background level to match the signal at the minimum at 
A(f) ~ 1. Center panel: three-particle correlation with the true background subtracted (corresponding to 
solid horizontal line in the left panel). Right panel: three-particle correlation with two-particle background 
normalized to the signal at A(f) ~ 1 (corresponding to the dashed horizontal line in the left panel). The 
maxima of the three-particle correlation plots in the center and right panels are set to the same value 0.1 in 
order to compare the details of the minimum regions. The three-particle correlation functions themselves 
have larger maxima, about 0.50 and 0.35 for the center and right panels, respectively. 

For illustration, we define the genuine three-particle correlation function to be the simple prod- 
uct of the two two-particle jet-correlation functions (i.e., the three-particle correlation function is 
factorized) : 

y3(A0i, A02) = y2(A0i)F2(A02). (40) 

This is shown in the center panel of Fig. [H In order to see the details of the minimum valley region, 
the plot maximum is set to 0.1; the three-particle correlation function itself has a larger maximum, 
about 0.50, which is on the near-side. Although in this specific, simple-minded example, the three- 
particle correlation cannot be zero if the two-particle correlation is non-zero, there are regions in 
the three-particle correlation function where the signal is almost zero as can been seen in the center 
panel of Fig. [H The two-particle correlation function, however, is significantly positive everywhere 
because one dimension is integrated out. Therefore, three-particle correlation is more powerful to 
determine the background than two-particle correlation. In principle this can be utilized to better 
normalize the background in two-particle jet-correlation. 

The three-particle signal after subtraction of the normalized (hence overestimated) background 
level is shown in the right panel of Fig. [H Again, in order to see the details of the minimum valley 
region, the plot maximum is set to 0.1; the three-particle correlation function itself has a larger 
maximum, about 0.35, which is on the near-side. As seen from the figure, the zero signal regions 
are now larger. For three-particle correlations more complex than the simple example here, the 
three-particle correlation function can have negative regions when the background normalization is 
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done with ZYAl or ZYAM for two-particle correlation. Again this is simply because one dimension 
of three-particle correlation is integrated out in two-particle correlation function. 

On the other hand, one may use three-particle ZYAM to obtain the background magnitude, i.e., 
the factor a, as aforementioned. This is a better estimation of the background because three- 
particle correlation contains richer, more differential information than two-particle correlation. The 
systematic uncertainty on a may be indicated by the change in the value of a from two-particle 
ZYAM to three-particle ZYAM. 

As discussed in Section 1X51 the value of h is determined from the associated multiplicity distribution 
in triggered events. The systematic uncertainty on h may be assessed by varying the size of the 
azimuthal regions used for three-particle ZYAM normalization. 



5 Summary 



In this paper, we have described an analysis method of jet-like three-particle azimuthal correla- 
tion between a high trigger particle and two softer particles. The method is devised with the 
particular two-component model in mind, in which the triggered event is made up by a di-jet and 
a bulk medium. The three-particle correlation analysis is extremely challenging due to the large 
combinatorial backgrounds in central relativistic heavy-ion collisions. Consequently, our emphasis 
is on the careful construction of the various backgrounds in an envisioned real data analysis and the 
evaluation of their major systematic uncertainties. By doing so, we demonstrate that three-particle 
jet-correlation result from such an analysis method is rather robust even though the background 
is large. The key to achieve the robustness is to use real data to construct the background terms 
self-consistently with only two well-defined scaling factors: the factor a to scale the inclusive event 
multiplicity to the background multiplicity underlying the triggered events by the ZYAl or ZYAM 
assumption, and the factor h to take into account any difference between the multiplicity distribu- 
tions of the inclusive events and the underlying background events such that the final three-particle 
correlation has zero yield at minimum. 
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